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Abstract 

We discuss constraints from cosmology on the models with direct gauge mediation based 
on the ISS supersymmetry breaking models. It is a generic feature that the ISS models possess 
a U(l) symmetry which is eventually spontaneously broken, when the models are based on the 
gauge symmetries with complex representations. We show that the resultant pseudo Nambu- 
Goldstone boson causes problems in cosmology, which leads to severe constraints on the 
models. We also show that most parameter space of the models can be probed by observing 
gravitational waves emitted from the cosmic string network when the U(l) symmetry is a 
gauge symmetry. 



1 Introduction 



The gauge mediated supersymmetry breaking models [U El El lU El El E] are one of the most 
attractive realizations of the phenomenologically viable minimal supersymmetric standard 
model (MSSM). So far, a lot of models of gauge mediation have been constructed, for the aim 
to improve the models in view of naturalness or consistency to cosmological observations. We 
have now many successful models to be examined by experiments. 

A downside of models of gauge mediations is, however, their complexity. The models 
are inevitably complicated, since we need to connect the supersymmetry breaking sector 
and the MSSM sector via gauge interaction of the MSSM. In this respect, models of direct 
gauge mediation where the MSSM gauge symmetries are embedded into flavor symmetries of 
the supersymmetry breaking sector are particularly attractive since the models require fewer 
particles and parameters (for a precise definition of direct mediation models, see Refs [8]). 
Before the discovery of dynamical supersymmetry breaking in simple supersymmetric gauge 
theories by Ref. [9], the construction of direct gauge mediation models could rely mostly on 
the vector-like dynamical supersymmetry breaking model developed in Refs. (HUH], the IYIT 
modelj^] Since it was shown that vector-like dynamical supersymmetry breaking is possible in 
more generic models by Ref. [9], many direct gauge mediation models have been constructed 
based on the ISS supersymmetry breaking models (See Refs. [131 H2 HSJ [El [T7] for early 
applications.)]^] 

There are two common features shared by the direct gauge mediation models based on the 
ISS supersymmetry breaking models. One, which is shared by all the direct gauge mediation 
models, is the Landau pole problem. Since we embed the MSSM gauge groups into flavor 
groups of the dynamical sector, there appear many fields which transform under the MSSM 
gauge groups. These fields contribute to the beta functions of the gauge coupling constants 
and drive them to Landau poles below the unification scale. This problem forces us to take 
the messenger scale rather high, while the supersymmetry breaking scale is kept rather low 
to realize the MSSM superparticle masses below the TeV range. The other feature shared 
is spontaneous symmetry breaking of global symmetries at the messenger scale. Especially, 
when the ISS models are based on gauge symmetries which have complex representations 
such as SU(N C ), the models possess a £7(1) global symmetry which is consistent with the 
mass terms of the models. We call it U(1)b symmetry in the rest of the paper. 

1 See Ref. for the successful direct mediation model based on this class of the vector-like dynamical 
symmetry breaking models. 

2 Another advantage of this class of models is that the models have simple realizations in string theory as 
emphasized in Refs. QH [H HOI HI] (see (22123] for reviews). 
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In this paper, we discuss the implications of these common properties, high messenger 
scale and spontaneous breaking of U{1)b symmetry, of the models. First, if the U(l)g is a 
global symmetry, we find that the pseudo Nambu-Goldstone boson (PNGB) resulting from its 
spontaneous breaking has a lifetime longer than the age of the universe. In this case, the energy 
density of the coherent oscillation of the PNGB easily dominates over the energy density of 
the universe. As a result, the models are consistent only when the reheating temperature 
of the universe is much less than hundreds GeV. Then, if the U{1)b symmetry is gauged, 
its spontaneous breaking does not predict the presence of Nambu-Goldstone boson but does 
generate cosmic strings with a tension of the order of the messenger scale. As we will discuss, 
most of the parameter space can be probed by observing gravitational waves emitted from the 
cosmic string networks in future experiments. Since the existence of the U(1)b symmetry is 
strongly interrelated to the structure of the gauge dynamics of the supersymmetry breaking 
sector, the observation of the gravitational wave gives us an important clue on the structure 
of the hidden sector. 

The organization of the paper is as follows. In section 2, we discuss the models of direct 
gauge mediation and spontaneous breaking of the global symmetries. In section 3, we discuss 
cosmological implications of the pseudo Nambu-Goldstone bosons resulting from spontaneous 
U(1)b symmetry breaking. In section 4, we discuss cosmological implications for the case 
where U(1)b symmetry is gauged. The final section is devoted to our conclusions. 

2 Direct gauge mediation 

In this section, we briefly review the models of direct gauge mediation based on the ISS 
supersymmetry breaking model. The Standard Model gauge groups are embedded into the 
global symmetries of the supersymmetry breaking sector. This class of models often involves 
global symmetries in addition to those identified with the Standard Model gauge groups. As 
we will discuss shortly, some of those extra symmetries are spontaneously broken at the scale 
of gauge mediation. The implications of the spontaneous breaking of those symmetries will 
be discussed in the following sections. 

2.1 An explicit example 

To be concrete, we discuss generic properties of the models of direct gauge mediation by 
considering an explicit model developed in Ref. [15], KOO model, as an example. The KOO 
model is based on the original ISS model [9], which consists of an SU(N C ) gauge theory with 
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Nf flavors, Qi and Qi, with a superpotential, 

W^Y^{m Q ) ij Q i Q j , (2.1) 

where mQ is a mass matrix. In the K00 model, the mass matrix rriQ is assumed to be 
diagonal but split into the first Nf — N c components with the eigenvalues being mo and the 
remaining N c components with /i , i.e. mq = diag(m , • • • , m , fio, ■ ■ ■ , fio). We assume that 
m o ^ A*o- I n this case, we may rewrite the superpotential into, 

w mass = m (g / g/) + ^{Q a Q a ) , (2.2) 

where J = 1, • • • ,Nj — N c and a = 1, • • • , iV c denote the flavor indices, while the color SU(N C ) 
indices are contracted in (QQ). 

In this model, the maximal global symmetry SU(Nf) is broken into SU (N = iV/ — iV c ) x 
SU(N C ) x Z7(l)p/ by the mass matrix. Here, U(l)p> denotes a £7(1) subgroup of SU(N C ) with 
the ratio of the charges of Q 1 and Q a being A^ c to N. In addition to those global symmetries, 
the model also possesses a U (1) symmetry with a charge assignment Q(+l) and Q(— 1), which 
we name [7(1)b, and a Z7(l) -R-symmetry with a charge assignment Q(+l) and Q(+l). 

To realize the gauge mediation mechanism, we embed the MSSM gauge groups into the 
subgroup of SU (N) x SU (N c ) global symmetries. By this, the MSSM sector is connected to 
the dynamical supersymmetry breaking sector in a direct manner. To construct a successful 
model, however, we need to break the .R-symmetry which forbids the gaugino mass terms. 
What is even more troublesome in the models based on the ISS models is that the gauginos 
do not obtain masses comparable to those of the sfermions even if the U(1)r symmetry is 
spontaneously broken as it happens in the model in Ref. [21], and hence, we need to break 
the R-symmetries explicitly^] 

In the KOO model, the comparable gaugino masses to the sfermion masses are achieved 
by adding the following term [15] , 

W dei = — -(Q'Qa^Qi) , (2.3) 
m x 

to the superpotential. Here, mx denotes a dimensionful parameter. With this additional 
term, the R symmetry is explicitly broken. The remaining global symmetries are SU(N) x 
SU(N C ) x U(l)p x [/(l) B n In the following discussion, we assume that mx is larger than 



3 For generic discussions on the gaugino masses, see Refs. [351 US]- See also Refs. [37J [551 HS1 GUI [HI] for 
related discussions. 

4 We define U(l)p as a linear combination of U(1)b and U(l)pi that has a trivial action on Qj. 
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m but much smaller than the Planck scale Mp^. Such higher dimensional terms can be, 
for example, generated by integrating out extra massive fields with masses m x coupling to 
(Qa, Qi) in a renormalizable theory. We also assume that other dimension four operators 
such as (Q a Qb)(Q a Q b ) are suppressed more significantly than those suppressed by m x , i.e., 
they are suppressed by the Planck scalej^] otherwise, the metastability of the vacuum could 
be spoiled. 

To explore the vacuum structure, it is useful to consider the magnetic description, where 
the model is well described by the composite fields 

Y I J = Q I Qj, Z I a = Q I Q a , Z a I = Q a Q I , <$>\ = Q a Q b . (2.4) 

In this description, the superpotential is given by 



W = hTi 



m 2 Y + /z 2 $ - xYx ~ X Z P - pZx - P^P - m z ZZ , (2.5) 



where p and x are components of magnetic quarks charged under the dual SU(Nf — N c ) 
gauge group. Here, we have suppressed all the indices of the gauge and global symmetries 
and defined 

m 2 = moA , p 2 = poA , m z = A 2 /mx ■ (2-6) 

Throughout this paper, we assume a dimensionful parameteij^] in magnetic theory to be the 
same order as A, i.e. h = A/A = 0(1). By appropriate phase rotations, we may take 
parameters m, p and m z to be real numbers without loss of generality. 



As long as the deformation superpotential (2.3) is small, the meta-stable vacuum identified 



in the ISS model still exists. In addition to the original supersymmetry breaking vacuum, 
this deformation also leads to other meta-stable and supersymmetry breaking vacua which 
are away from the ISS vacuum. These vacua are given by 

n = ^fe);+o(m z ), 

X 1 j = mS*, , 





= — (nr 


m z NK 


x T j 


= mSj , 


P\ 


= P Tl a , 




mp r i 




m z a 



J ' 



P a i = PT' 



(2.7) 



i j 



~Z\ = -^T% 

m, 



5 Such suppressions can be explained if we assume an approximate i?-symmetry with the charge assignment 
R{Qi) — R{Qi) — 1 ancl R{Qa) = R(Q a ) — 0, which is explicitly broken only by the fJ,Q a Q a term [15] . 

6 The dimensionful parameter A is related to the scale in electric description A and the one in magnetic 
description A by 

^3N C -N f ^3{N f -N a )-N f = f^Nf-N^Nf _ 
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where n is any integer between and N. The matrices T a j and Y I a have 1 in the first n 
diagonal elements and elsewhere; 

pa _ j O nx (Ar_ n ) \ _ / E„ O nx (Ar c _ n ) \ ^ 

' \0(Ar c -n)xn ®(N c -n)x{N-n) ) ' " \0(N-n)xn 0(N-n)x(N c -n) J 

The matrix l m is an m x m identity matrix and 'SF is a g x q matrix whose first p diagonal 
elements are 1 and otherwise. The vacuum with n = corresponds to the vacuum of the 
ISS model, while the vacua with n > have lower energy than the ISS vacuum, and the 
supersymmetry breaking scale of these vacua is 0(p 2 ). As we show in the appendix A, all 
pseudo-moduli are stabilized by the Coleman- Weinberg potential. 



At the vacuum in Eq. (2.7), the global symmetries are spontaneously broken into SU(n) x 
SU(N — n)x SU (N c — n) x U (l) 2 for n > 0. Here, two preserved {7(1) symmetries are linear 
combinations of the following three U(l) symmetries; 

U(1) P : p -> e Wp p . (2.9) 
U(l) c :p^( el9C l nXn i^ Jn ° ) P , (2.10) 

\ U e N o-n l (A r c _ n)x(A r c _„,) J 

U{l) d :p^ P \ iapfc^,. ° I • (2-11) 

Then one of the preserved symmetries is a combination of U(l)p and U(l) c with Op = —9 C 
and the other is a combination of U(l) c and U(l)d with 9 C = —9d- In contrast, at the highest 
energy vacuum, i.e. the ISS vacuum, only a U(l) global symmetry is spontaneously broken 
and the resulting preserved symmetry is SU(Nf — N c ) x SU(N C ) x U(l)p. In both cases, one 
of the i7(l) symmetries is spontaneously broken, which corresponds to the U(l)p symmetry 
in the electric description. 

Finally, let us embed the MSSM gauge groups into the flavor symmetry and construct 
direct gauge mediation. In this case, the admixtures of the dual quarks and the mesons play 
roles of the messengers whose masses are around m. The detailed analysis of the mass matrix 
of the supersymmetry breaking sector including the messengers is given in the appendix 
A. When we embed the MSSM gauge groups into SU(N — n), the messenger mass matrix 



squared of the fermions is given by A 2 in (A.9). This is nothing but the one appearing in 



O'Raifeartaigh model of KOO vacua [15]. Thus, just replacing N c , N in the result of [15] with 
N c — n, N — n respectively, we obtain the leading order of the gaugino masses, 

q 2 (N c — n)^ (9, , , , q 2 (N c — n) ha 2 m z 
m x ~ \ A \ 2 ' F*— logdetM = y \ c ' £ . 2.12 

(47r) 2 9$ (47r) 2 m z ($)-m 2 
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Here, (<£>) is determined by the Coleman- Weinberg potential for full O'Raifeartaigh model 
including the contributions from A\ and B\ (see the appendix A), which is proportional to 
the explicit R-breaking parameter m z , i.e. ($) ~ m z . As a result, at the leading order of 
m z /m, the gaugino masses are given by, 

When we embed the MSSM gauge groups into the SU (N c — n) group, gaugino masses are 
given by 

Notice that if we embed the MSSM gauge groups into the SU(n) symmetry, the messenger 
sector does not include a tachyonic direction anywhere in the moduli space, and hence, one 
expects that gaugino masses vanish [25] • in fact, we can show detM = fi 2 m z by a direct 
computation, which indicates that the determinant of the matrix is independent of $, and 
hence, gaugino masses vanish. 

The sfermion masses are, on the other hand, roughly given byQ 

9 2 ^ fn 1 K\ 

m sfermion ~ ~. • K 1 - 1 ^) 

All m 

In order to have comparable gaugino masses to the sfermion masses in the TeV range, the 
mass parameters should satisfy 

u 2 

~ 100 TeV, (2.16) 

m 

m z ~ m . (2-17) 
Here, the second condition is realized when 

A ~ y/mm^ . (2.18) 



As discussed in Ref . [To] , the Landau pole problem is solved by assuming that the messen- 
ger scale is higher than 10 12 GeV. Besides, the supersymmetry breaking scale is restricted to 
be smaller than 10 9 ' 5 GeV to avoid unacceptably large flavor violating soft parameters from 



7 The sfermion squared mass contributions from x's correspond to the ones discussed in Ref. [521 1551 151) . 
which tend to be suppressed compared with the other contributions from p's and Z's. The detailed analysis 
on the sfermion squared masses will be given elsewhere. 
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the gravity mediation effects. By combining these constraints with Eq. (2.16), we find that 
the model is viable for 



10 12 GeV < m < 10 14 GeV , 

10 8 - 5 GeV< fi <10 9 ' 5 GeV. (2.19) 

Let us summarize the important features of the models with direct gauge mediation based 
on the ISS supersymmetry breaking models. In the KOO model, the global U(l) b symmetry is 



spontaneously broken by the expectation values of x an d X i n Eq. (2.7). Since the messengers 
obtain masses from the vacuum expectation value of x an d Xi the U{1)b breaking scale 
coincides with the messenger scale. Notice that this non-trivial interrelation between the 
messenger scale and the U(1)b breaking scale is rather generic feature of the models of direct 
gauge mediation based on the ISS variant models. 

In the following discussion, we concentrate on the vacuum with n = 0, since it is the 
vacuum closest to the symmetry enhancement point which can be dynamically chosen in 
the cosmological evolution]^] We also assume the models where the MSSM gauge groups are 
embedded into SU (N) with an assumption iV c > 10 to avoid the Landau pole problem (see 
[15] for details). 



3 Cosmological implication of the PNGB 

As emphasized in the previous section, spontaneous symmetry breaking of global symmetries 
is expected at the messenger scale in the direct gauge mediation models. If the symmetries 
are exact, the corresponding Nambu-Goldstone bosons are exactly massless. It is generally 
believed, however, that any global symmetries should be broken in gravitational theoryj^] That 
is, the global symmetries are expected to be explicitly broken at least by higher dimensional 
operators which are suppressed by the Planck scale. Therefore the resultant Nambu-Goldstone 
bosons are pseudo Nambu-Goldstone bosons (PNGBs) and obtain finite masses. 

When the PNGBs are massive but very light, their lifetime can be very long and cause 
problems in the cosmological evolution. To avoid such problems, the symmetry breaking 
operators should not be suppressed so much so that the resultant PNGBs are heavy enough. 
As we will see, however, the U(1)b symmetry is highly protected by the gauge symmetry 
and the corresponding PNGB cannot be heavy enough] 1 " 9 ] As a result, the PNGB from U(1)b 

8 We may straightforwardly extend our analysis to those with n > 0. 

9 This can also be understood by the fact that there is no global symmetry in string theory. See |35j for a 
recent discussion. 

10 At the vacuum with n > 0, we also have spontaneous breaking of the non-abelian global symmetries. 
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breaking dominates over the energy density of universe, which ruins the success of the standard 
cosmology. 



3.1 PNGB mass 

The distinctive feature of the U(1)b symmetry is that it is highly protected by the SU(N C ) 
gauge symmetry of the supersymmetry breaking dynamics. Even the lowest dimensional 
operators which explicitly break the U(1)b symmetry are quite suppressed; 



w (leading) C B n N c C§ ^ Nc ,, 

+ J^ Q ' (3 - 1} 



where Cb b are dimensionless coefficients of 0(1). These operators correspond to the baryon 
and anti-baryon operators, respectively. By remembering that the size of the gauge group 
SU(N C ) is required to be large, N c > 10 (see Ref. [TS] for details), we see that these operators 
are quite suppressed. 

In the magnetic dual description, the above baryon and anti-baryon operators are trans- 
lated into, 

\2Nc-N* \2Nc-Nf 

~ ^X"^ + eloper*"'-"- • (3-2) 

lvlp L 1V1 PL 

where c Bj § are again the dimensionless coefficients of 0(1). The dynamical scale A of the 
SU(N C ) gauge theory has been incorporated based on the dimensional analysis. 



By combining the superpotential of the K00 model in Eq. (2.5) with the above U(1)b 
breaking operators, we obtain, 

A 2 ~ A 2 

W = mAM + qqM + ZZ + - M 2 + W B + -- - , 

m x M PL 



r , - ^(leading) M A^^X^M A ^-N f+ l~N f -N c . M 

"u + ^ + J^l +•--, (3-3) 



where the meson field M collectively denotes the meson fields Y, $ and Z, and the dual 
quarks, q and q denote q = (x, p) and q = (x,p). We have omitted detailed flavor structures 
of each term. In the above superpotential, we included the higher dimensional U(1)b breaking 
operators which are relevant for the PNGB mass. 

The leading contribution to the PNGB mass is from the cross term of the fourth term of 



W and the second and third terms of Wb in Eq. (3.3) via C ~ I^mW 7 ! 2 , which roughly leads 



Those symmetries, however, can be broken even by the tree-level superpotential in Eq. (2.1 1. Therefore, the 
PNGBs from the non-abelian global symmetries may have masses at the messenger scale. 
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Figure 1: Left) The parameter region of the dynamical scale A (see Eq. (3.9)). Right) The PNGB 
mass as a function of m. The upper solid line corresponds to ttt-pngb for A = A max for N c = 10, 



where ttt-pngb is dominated by the one in Eq. (3.4). The lower solid line corresponds to ttt-pngb f° r 
for A = A m i n for N c = 10, where mpNGB is dominated by the one in Eq. (3.7). The dashed lines 
denote those values for N r = 11. 



to, 



^2N c -N f +3 



m 



PNGB 



M 



PL 



m N f ~N c -l 



(3.4) 



Here, we have used (<&) ~ m z ~ m. Notice that the contributions from the second and the 
third terms of W are vanishing at the leading order since (Zx) = (xZ) = 0. 



Another relevant contribution comes from the A-terms of W 



(leading) 
B 



sional Kahler potential term, 



5K 



MM^qqi 
A2 



via the higher dimen- 



(3.5) 



By plugging M = m + 9 fx and renormalizing the kinetic term of q, we obtain the A-term of 

(leading) 



which leads to 



C~( Nf -N c )^mxW B 



(leading) 



+ h.c. 



l PNGB 



(N f - N c )li- 



2N c -N f -2 



M 



N c -3 
PL 



-m 



N f -N c -1 



(3.6) 



(3.7) 



By comparing Eqs. (3.4) and (3.7) with the conditions in Eq. (2.16), we find that the PNGB 
mass from the A-terms is dominant when the dynamical scale L is relatively low (see Figurejl] 
for numerical results). 
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Before giving the numerical results of the PNGB mass, let us discuss the constraints on the 
dynamical scale A. For generating the gaugino masses comparable to the sfermion masses, the 
.R-breaking mass term, m^ZZ, played an important role. As we mentioned earlier, the other 
mass terms such as $ 2 should be, on the other hand, much smaller than m^ZZ otherwise the 
metastable supersymmetry breaking vacua would be destabilized. Thus, we need to require 
mx <C Mpl- By remembering that we also require A ~ y/mmx for a successful gauge 



mediation (see Eq. (2.18)), this condition leads to the upper bound on the dynamical scale 
for a given messenger scale. 

The dynamical scale A is also constrained from below, A ^> m, otherwise the metastability 
of the supersymmetry breaking vacua could not be analyzed. That is, the incalculable higher 
dimensional terms of the Kahler potential of the pseudo-flat direction such as, 



1$ 



4 



K ~ ^ , (3.8) 

become comparable to the radiatively generated Kahler potentials via the dual quark inter- 
actions. 

In Fig.[TJ we show the upper and lower limits on A as a function of m which are defined 

by, 



A max = y/0.1 x M PL m , A min = 10 x m . (3.9) 

The figure shows that the dynamical scale is in a range of 10 13 ~ 16 GeV. In the right panel 
of Fig.[l| we show the PNGB mass for a given messenger scale. The figure shows that the 
PNGB mass is quite suppressed, and for example, 

m PNGB <0(l)keV, (3.10) 

for m ~ 10 13 GeV. As we will see shortly, the energy density of the PNGB dominates over 
the energy density of the universe due to the long lifetime, although it is very light. 

Before closing this section, let us discuss the masses of the fermionic and the scalar super- 
partners of the PNGB. Since the U(1)b symmetry is broken at the energy scale higher than 
the supersymmetry breaking scale, these partners are also expected to be light. Unlike the 
PNGB, however, they obtain masses by the supersymmetry breaking effects. To discuss their 
masses, it is useful to define the PNGB chiral supermultiplet A by, 



A -A 



X ~ m x , x~mxe m . (3.11) 

The partners of the PNGB obtain the masses via the Kahler potential, 

h 2 $t$ 

K~^- 2 —(A + Ar + -- - , (3.12) 
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which is radiatively generated by the \ ~ % — p interactions. By plugging $ ~ m + p 2 8 2 into 
the above Kahler potential, we obtain the masses of the superpartners of the PNGB; 

h 2 P 2 h p 2 

m sca l ar ~ . (3.13) 



An 2 m An m 

Therefore, they are typically heavier than the superparticles in the MSSM unlike the PNGB 
of the i7(l) b symmetry. 

3.2 PNGB decay width 

As pointed out in Ref. [36J, discrete symmetries such as charge conjugation symmetry and 
CP-symmetry play important roles to determine the interaction between the PNGB and the 
MSSM sector. 

Let us first consider a charge conjugation transformation defined by the exchange, 

Q^Q , (3.14) 

with appropriate charge conjugations of the gauge fields. The PNGB is odd under this 
transformation. This charge conjugation is a good symmetry when we neglect the matter 
fields in the MSSM sector. Thus, as long as we integrate out only the messenger fields, the 
effective interactions of the PNGB respect this symmetry. As a result, the lowest dimensional 
operators which are relevant for the decay of the PNGB into the MSSM gauge fields require 
at least three gauge fields which are roughly given by [371 EE] , 

£cfr ~ 9 ^^{D p F aP ){D^F aT ){D"D a F^)V , (3.15) 
Air m' 

where V denotes the PNGB, F pu the gauge field strengths of the Standard Model, and _D M the 
appropriate covariant derivatives. We collectively denoted the corresponding gauge coupling 
constants by gsM- Therefore, the decay width of the PNGB into photons through this operator 
is negligibly small. 

Once we integrate the MSSM sector as well, the above non-renormalizable operator could 
lead to the effective interactions between the PNGB and the Standard Model fermions, /, at 
the three-loop level, which are at most given by, 

£cff ~©^ WV/ + ' ,,) ' (3 ' 16) 

Here, the gauge coupling constants are estimated at the messenger scale. The decay widths 
of the PNGB into the fermions (i.e. the electrons or neutrinos) through these operators are 
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quite suppressed, and hence, the dominant process should be the one into photons through 
the diagrams with one more loop. As a result, the decay width is at most of the order of 

V 47r / m 2 V IkeV J \ m J 

Therefore, we find that the lifetime of the PNGB is much longer than the age of the universe p] 
The PNGB may also decay via the interaction terms which explicitly break the U(1)b 



symmetry in Eq. (3.3). From the explicit breaking term W^ eadmg \ the messenger fields obtain 



an effective mass term which depends on the PNGB, 

A 2N C -N f 



6m q ~ C 6 Wr^-V^ . (3.18) 

M PL 

Thus, after integrating out messengers, we obtain an interaction terms, 

1 r 1 K2Nc-N f 

— 2 / d6* log M me88 W a W« ~ Tg-jO* - c^—^m^-^VF^ . (3.19) 

Notice that these operators vanish for c# = c§ since the charge conjugation is a good sym- 
metry in this limit. As a result, the PNGB decays into two photons via these operators with 
the decay width, 



1 m PNGB 

16tt M e 2 ff 



(3.20) 



where we have defined 



Mgs = -i^M Pi (^- N ' (*™) W » Mpl . (3.21) 
c B -c B V A / \ m J 

The last inequality is obtained by remembering Nf — N c > 5 and 2N C — Nf > 0. As a result, 
the decay of the PNGB via the U(1) B breaking operators is subdominant, and hence, the 
lifetime of the PNGB is much longer than the age of the universe. 

3.3 PNGB dominated universe 

Now we are ready to study energy density of a coherent oscillation of the PNGB in the 
history of cosmology. During inflation, the U(1)b symmetry is expected to be restored by the 
inflation dynamics if the Hubble scale during inflation is larger than m, since x's are expected 
to obtain effective masses around the symmetry enhancement point, i.e. x = X = 0- I n this 



L The scalar and the fermionic partners of the PNGB decay immediately into a pair of the PNGBs or a 



pair of the PNGB and the gravitino respectively via the operator in Eq. (3.12) 
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case, the initial amplitude of the coherent oscillation of the PNGB from its true minims^] 
after inflation is 0(m). If the Hubble scale during inflation is smaller than m, the U(1)b 
symmetry is not restored during inflation. Even in this case, the initial amplitude of the 
coherent oscillation of the PNGB after inflation is again expected to be 0(m), since the 
PNGB has a very flat potential on which the initial condition is randomly distributed. Thus, 
as long as the Hubble scale during inflation is much larger than wpngbi the initial amplitude 
of the coherent oscillation of the PNGB after inflation is expected to be 0(m). Thus, the 
initial energy density of the oscillation of the PNGB after inflation is expected to be, 

Ppngb ~ mp NGB m 2 . (3.22) 

In the followings, we discuss the evolution of the energy density of the PNGB after inflation 
with this initial condition. 

First, let us assume that the reheating temperature of the universe after inflation is lower 
than 

^osc ~ V m PNGBM PL . (3.23) 

In this case, the PNGB starts oscillating around its minimum when the Hubble parameter 
becomes lower than 

H < H osc ~ mp N GB , (3.24) 

which is before the radiation dominated era (see Fig.|2|. Before the universe enters into the 
radiation dominated era, both the inflaton and the PNGB oscillate, and hence, their energy 
densities scales in the same way by the cosmic expansion. Thus, at the beginning of the 
radiation dominated era, the energy density of the PNGB at the reheating time is given by, 

PPNGB — PPNGB X ~7c) — , , 2 • (3.25) 
P} 1V1 PL 

Here, pi denotes the energy density of the inflaton, Tr the reheating temperature, and the 
superscripts of the energy densities correspond to the ones in Fig.[2j In the final expression, 



we have used Eq. (3.22) and 



pf ] ~ m PNGB MpL , 

P? ) ~ n, (3.26) 
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Let us remember that the U(1)b symmetry is explicitly broken. 
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log V 



PpNGB 






A* 




PSM 




Tosc T R Tdom -logT 

Figure 2: The evolution of the energy densities as the universe cools down. The two solid lines 
represent the density for inhatons and that for pseudo-Nambu-Goldstone bosons, respectively, and 
the dotted line represents that for Standard Model matters. To avoid spoiling the successful Standard 
Cosmology, the energy density, ppngb, should not exceed psm at the temperature above T eq = 
0(1) eV (see text). 



where the first equality is obtained from Eq. (3.24) and the second one represents the energy 
density conservation at the reheating time. 

During the radiation dominated era, the energy density of the radiation scales by the 
fourth power of the inverse of the scalar factor a, while the PNGB energy density scales by 
the third power of the inverse of a. Thus, by remembering that the cosmic temperature scales 
the inverse of a, we find that the PNGB energy density dominates the energy density when 
the temperature becomes lower than, 

(B) m 2 

Tdom ~T R x — ~T fi x . (3.27) 

P R M PL 

Since the lifetime of the PNGB is longer than the age of the universe, the domination tem- 
perature T dom should be much lower than the temperature of the matter-radiation equality 
temperature of the Standard Cosmology, T oq = 0(1) eV. Thus, for the successful cosmology, 
the reheating temperature is severely restricted to 

T R ~ T dom x « Tcq x ¥k „ 0( ioo) GeV x f 1 !^!) ' . (3 . 28) 

This constraint is rather severe when we consider baryogenesis mechanisms which require very 
high reheating temperature such as the thermal leptogenesis mechanism [39J . 

We may do a similar analysis for T osc <T R . In this case, the PNGB oscillation dominates 
over the energy density when the cosmic temperature becomes lower than 

2 

Tfl 

T dom ~ T osc x . (3.29) 

1V1 PL 
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Thus, by requiring T^ om <C T eq , we find that 



/10 13 GeV\ 4 

mpNGB < 10- 10 keV x . (3.30) 

\ m J 

By comparing the PNGB mass in Fig.[TJ we find that this condition is difficult to be satisfied 
without causing the Landau pole problem. 

4 Gauged U(1)b symmetry 

A simple way out of the above discussion is the gauged U{1)b symmetry. When the U{1)b 
symmetry is a gauge symmetry, its spontaneous breaking does not lead to the Nambu- 
Goldstone boson, and the cosmological constraints discussed above is not applicable. Instead, 
it leads to another cosmologically interesting object, the stable cosmic string. In this section, 
we discuss the implications of the cosmic string in the direct gauge mediation models based 
on the deformed ISS models. 

4.1 Classification of vortex solutions 

Here, we show that there are two types of vortices in the deformed ISS-model on which the 



KOO gauge mediation model bases ^ The solutions corresponding to the topological defects 
have finite energy, and hence, the configuration of the solutions decay rapidly enough as 
the spatial radius increases. As is well known, this condition is satisfied when the scalar 
fields responsible for spontaneous symmetry breaking take values on the same gauge orbit at 
infinity (see Ref. [H] for a review.) For example, the one dimensional topological defect, i.e. 
the vortex, corresponds to a map from the spatial infinity of two-dimensional space to the 
gauge orbit which defines the winding number 7ri(Gi cai/-ffiocai) (see also the Appendix B). 
The homotopy group, vr 1 (Gi oca i/ifi oca i), can be calculated by using the exact sequence, 

tt 2 (G) = 4 n 2 (G/H) 4 TTi(ff) A 711(G) 4 71*1 (G/ H) 4 7t (H) = , 

Ker/fc+i = lmf k . (4.1) 

Let us start with the highest energy vacuum where the only non- vanishing dual quarks in 



(2.5) are x 1 j = X 1 j = m & j U P to symmetry transformations. In this case, symmetry breaking 



occurs as follows: 

G = SU(N) color x SU(N) x SU(N C ) x U(l) B x U(l) P ^H = SU(N) diag x SU(N C ) x [7(l)j 



13 Topological defects in the ISS-metastable vacuum were firstly studied in Ref. [40] . where the semi-local 
vortex solution were explicitly constructed. 



15 



Since U(l)p acts trivially on this metastable vacuum, non-trivial contribution comes only 
from the gauged U(1)b breaking, 

« (^)=«( suw * um °) = z. 



H, 



local 



Therefore, we find that there is at least one finite energy vortex. 

The vortex solutions at the lower energy vacua are more interesting. Gauging only U(1)b 
is very similar to the case of the highest KOO vacuum. So we here consider a vortex when 
both U{1)b and U{l)p are gauged. The fundamental group is 

(mm - urns x u{i) r • z x z _ 



So there should be two types of finite vortices. As we discuss in the Appendix B, one of the 
vortices corresponds to the so-called semi-local vortex. In the following analysis, we again 
concentrate on the case of the highest energy vacuum, n — 0. 

4.2 Cosmic strings 

Now let us discuss the cosmological implications of the vortex solutions. Since the U(1)b 
symmetry is gauged and vr 1 (Gi oca i/ff loca i) = Z, dynamical generation of the cosmic strings 
with a tension 

/i T ~7r 2 m (4.2) 

is expected at the phase transition of the U(1)b symmetry. Such phase transition occurs 
when the Hubble scale or the temperature of the universe become lower than the symmetry 
breaking scale m. At the phase transition, typically one cosmic string is generated in a Hubble 
volume 



ii 



First, let us consider the case with m > if in f,Tff, where H [n f is the Hubble scale during 
inflation and Th the highest temperature afte inflation 15 In this case, the cosmic strings 
are diluted away by the inflation even if they are generated at the phase transition before 
the inflation. Therefore, we find that the models with direct gauge mediation based on the 
deformed ISS is consistent with cosmology if the U(1)b symmetry is a gauge symmetry, since 



14 This estimation represents that the causality does not permit an exchange of information beyond the 
horizon scale. In Ref. [43], however, it is pointed out that if the system cools quickly a field value can be 
different within a causally connected region. As recently emphasized in Ref. [44] . this Kibble-Zurek mechanism 
provides a substantially larger abundance of topological defects from phase transitions in the early universe 
than the original estimation by Kibble. Notice however that the energy density of the cosmic strings on its 
simple scaling low (see below) is insensitive to the initial energy density of the cosmic strings. 

15 The highest temperature Th during inflation is given by Th — (Tj-ff; n f Mpl) 1 / 4 . 
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cosmic strings do not leave any visible effects after inflation. Since the messenger scale m is 
above 10 12 GeV, the condition m > Tn is still consistent with thermal leptogenesis [39J. 

Next, let us consider the cases with Max [if; n f,T#] > m where the cosmic strings are not 
inflated away. In these cases, the phase transition takes place after the end of inflation, and 
the produced cosmic strings survive until now. 

Once very long cosmic strings are formed at the phase transition, their energy density is 
naively expected to scale as p str ~ fiTH 2 a~ 2 , where is the Hubble parameter at the phase 
transition. If this is the case, the energy density of the cosmic strings eventually dominates 
over the energy density of the universe unless the tension fij< is negligibly small. However, 
the long cosmic strings break into loops through the reconnection processes of the strings. 
The produced loops subsequently disappear by emitting gravitational waves. In this way, the 
energy density of the string finally reaches the simple scaling, p str oc HtH 2 , despite the very 
complicated system of the string network (see jl5] for details). As a result, the energy density 
of the stable cosmic string is always subdominant at any time as long as the tension \xt is 
much smaller than the Planck scale. 

One of the striking effects of the cosmic string is its contributions to the anisotropy of the 
cosmic microwave background radiation (CMBR). Recent CMBR observations from WMAP 
and SDSS experiments have placed a limit on the tension of the cosmic string, G\xt < 3.5 x 
1CT 7 [46J where G is the Newton constant, G ~ (1.2 x 10 19 GeV) 2 . In terms of the messenger 
scales, the above constraint corresponds to 

m<10 15 - 2 GeV. (4.3) 

Therefore, the messenger scale we are interested in is consistent with the constraints from the 
anisotropy of the cosmic background. 

Another striking signature from the cosmic strings is the gravitational radiation produced 
by the string loops [17J H3] . The most stringent constraint on the cosmic strings through the 
gravitational wave observation is the one coming from the pulsar timing tests. According 
to Ref. [52], the tension of the cosmic string is constrained to be G^t < 10 -9 by the pulsar 
timing limitsj^] Thus, we find that the messenger scale of the direct gauge mediation is also 
consistent with the constraint obtained from the pulsar timing limits. 

16 This value is the result when we assume that the typical length of the newly produced loops at the Hubble 
scale H is about Aoop ~0.1x H^ 1 . This length at the production is suggested by numerical simulations done 
in [?5J HH1 \5U\ . If the typical length at the production is much shorter, the pulsar timing test is not sensitive 
to the gravitational wave from the loops of the cosmic strings. In such cases, the gravitational waves from 
the infinite strings become important at the frequency of the gravitational wave around 10 -8 Hz to which the 
pulsar timing test is sensitive [51] . 
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The future experiments of the gravitational wave observation will be much more sensitive 
to the gravitational waves from the cosmic strings. As discussed in Ref. [52], we may probe 
the gravitational waves from the cosmic strings with the tension down to GfiT — 10~ 16 by 
the sensitivity of the LISA experiment for example p] Thus, the future gravitational wave 
experiments cover most of the messenger scale, 

m > 10 10 - 5 GeV . (4.4) 

Therefore, if the phase transition of the U(1)b takes place after the end of inflation, the 
models based on the deformed ISS can be tested by observing the gravitational wave. 

Finally, let us comment that the Hubble scale during inflation can be determined experi- 
mentally by the tensor to scalar ratio r of the cosmic microwave background, since they are 
related by ifi n f ~ 10 14 x r 1//2 GeV. The COrE experiment |53j, for example, aims to measure 
the ratio down to 10~ 3 which corresponds to ifi n f ~ 10 12-13 GeV. In this way, we can indirectly 
test one of the conditions of the occurrence of the phase transition after the end of inflation, 
i.e. H [n { > m, via the observation of the tensor mode of the cosmic microwave background. 

5 Conclusions 

In this paper, we discussed cosmological implications of the models with direct gauge medi- 
ation based on the deformed ISS models. The generic properties of this class of models is 
that some of the global symmetries are spontaneously broken at the messenger scale. Es- 
pecially, we investigated the implications of spontaneous breaking of the U{1)b symmetry 
which is an intrinsic feature of the ISS models based on the gauge symmetries with complex 
representations. 

When the U(1)b symmetry is a global symmetry, it is expected to be broken explicitly 
by gravitational interactions. In the ISS models, however, the U(1)b symmetry is highly 
protected by gauge symmetry of the dynamics of supersymmetry breaking. Due to such 
protection, the mass of the resultant PNGB is very small and its lifetime is longer than 
the age of the universe. The energy density of the coherent oscillation of the PNGB easily 
dominates over the energy density of the universe. As a result, we found that the models 
with the global U(1)b symmetry are consistent only when the reheating temperature of the 
universe is much less than hundreds GeV. This constraint is rather severe when we discuss 
baryogenesis mechanisms such as leptogenesis which requires the reheating temperature higher 
than T R > 10 9 GeV [39]. 

17 This number is again for £\ oop ~ 0.1 x H^ 1 at the production. See Ref. [52] for details. 
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When the U(1)b symmetry is a gauge symmetry, on the other hand, its spontaneous 
breaking predicts the cosmic strings instead of the Nambu-Goldstone boson. We found that 
the messenger scale of our interest, i.e. m = O(10 12 ~ 14 ) GeV, is consistent with the cosmic 
microwave background constraints and the pulser timing limits even if the phase transition 
occurs after inflation. Furthermore, we also found that most of the parameter space can be 
tested by observing the gravitational waves from the cosmic string networks at the future 
experiments. 

Finally, let us emphasize again that the existence of the U(1)b symmetry is strongly 
interrelated to the structure of gauge dynamics of the supersymmetry breaking sector. Thus, 
the detectability of the remnant of spontaneous breaking of the U(1)b symmetry allows us 
to get a glimpse of the structure of the hidden supersymmetry breaking sector through the 
observation of the gravitational waves. The U(1)b symmetry is, however, not expected in 
the IYIT supersymmetry breaking models since they are based on the gauge theories with 
real representations. The supersymmetry breaking vacuum of the IYIT model is stable, while 
that of the ISS model is metastable. Therefore, by investigating the implications of the U(1)b 
symmetry breaking, we will not only be able to put constraints on the model parameters, but 
also might get a hint on a speculative question; whether we are on the metastable vacuum or 
not. 
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Appendix A Stability of lower vacua 

In this appendix, we study the stability of the lower energy uplifted vacua in the model studied 
in [15] by calculating the Coleman- Weinberg potential. The supersymmetry condition for the 
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superpotential (2.5) dW = becomes 



p 2 - pp = , 
™ 2 ~ XX = , 

Yx + Zp = J X Y + P Z = 0, (A.l) 
PX + m z Z = , xp + m z Z = 
Zx + $p = Q, xZ + p® = Q. 
Except the first equation, all equations are satisfied, and we obtain supersymmetry breaking 



vacua. Note that there are flat directions. This can be seen by noting that (A.l) is still 
satisfied by the change of fields 

X ->■ AhxBn , x ->■ B^xAn , 
P -> A^pC Ne , p -> Cj£pA N , 

_! ~ „_V ( A - 2 ) 

^ ~~ ^iv ^^Vc ' ^ — * C Nc ZB N , 
$ -> C£$CW C , F -> S^FSjv , 

where A, 5, C denote invertible matrices whose sizes are shown in subscripts. When A, B, C 
are especially special unitary matrices, they generate Goldstone modes. This explains half of 
the possible non-trivial modes. The other half are pseudo-moduli and the quantum fluctuation 
may develop a non-trivial Coleman- Weinberg potential. 

It suffices to consider the following background fields to obtain the Coleman- Weinberg 
potential. 





I, 2 


2 

foa m fun \a , / ^iN c -n-, a 

$ b = — (livj b + 7(1 N c ) b 


X*j 


= amSj , 


X J j = a^mdj , 


p\ 


= afipT^ , 


p a i = c^/rVn , 




(3m p j 
m z a 


~ a P~ l mp^ a 

> a i — 1 i > 

m z 



(A.3) 



where a, (3 and 7 are pseudo moduli. (1'^ ") b is an N c x N c matrix whose lower N c — n 
diagonal components are 1. The boson and fermion mass matrices are given by 



., :w^ c w ch w^ abc w c \ 2 _(w^ c w ch , 
Mi^ir ■ \\ m \v'j ■ m r-{ w ar w ]ch] - 

From this, the Coleman- Weinberg is computed by using the formula 



V eff = -L <Tr m% log ^| - Tr m% log . (A.5) 
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Components that are not directly coupled to W P p$W^ do not contribute to the Coleman- 
Weinberg potential due to the boson-fermion cancellation. Among the off-diagonal compo- 
nents W a j, c W' c of the boson mass matrix, the only non-vanishing component is 

Therefore we focus on the components of the boson mass matrix m 2 B that are connected 
to Wpp^W^. It is useful to divide the index a into a\ and a 2 , where a\ = l,-- ,n and 
a 2 = n + 1, • • • ,N C . Similarly, divide I so that / = (h,I 2 ) such that 1\ = l,-- ,n and 
I 2 — n + 1, • • • ,N. Then we can check that those that are coupled to the off-diagonal element 
Wpp$W^ can be written as 



where 



m Seffective 




J) 


® Hn(iV c - 


-n) £ 




B 2 \ 
A 2 J 


® ~R(N-n)(N c ~n 


)] , (A.7) 






$7 

J2 








Afi 
e 2 












Zj M 1 


\ 


01 




2 
















a 2 /3mfj, 


$t&i 
















a/3 









mfi 

a 2 /3 




p\T 















mm > + ar/m 











a/3 







a 2 / 


? 2 /z 2 + a 2 m 2 + 7 2 





amm z + — 

z a 






a 2 (3fim 





mniz 
a 


-f 07m 









a 2 m 2 + m 2 




Z\h 


V 




rrt/i 

a 2 /3 

re 2 









amm z + — 

^M 2 





(A.8) 


/ 




p\T 


/ 


# + 7 2 









mm > + orfm 





\ 






A 2 = 


_ « a 





a 2 m 2 + 7 2 





a 2 m 2 + m 2 


amm z + — 

z a 




(A.9) 








V 





amm z + 


1117 

a 





















ei 


2 $t/i 

e 2 




^e 2 


ry\e 2 
Z Mi 


e 2 














$t&2 
Oi 
















^ 














$t&i 
a 2 






























tl2 
_ P/l 

1 C12 












-P 2 













J 


(A.10) 










































012 


V 














) 









21 



Bo 





\ r 


ptM 2 
l J e 2 


7^2 
M 2 


y\M 2 

e 2 




Ph 


( 


-P 2 








\ 


P^ 















z\T 


















V o 











/ 



(A.11) 



Note that the first summand of the boson matrix matrix (A. 7) always has a zero eigenvalue 
whose two eigenvectors are 



m 2 — 7 A a/3 A „ m 2 — 7 A 

,0,^,0,-/3,0,0, *-,0, 



m 



rn/j, 



A 1 

apm p 



A 



m 2 — 7 A 

0, -,o, , 

apm 



,0, 



1 m 2 -7A o/3A 



rrifi apm p mp: 
where one is a Goldstone mode and the other a pseudo-modulus. 



-/3,0 



(A.12) 



Using the mass matrices and the formula (A. 5), we can compute the effective potential. 



Note that the mass matrix expressions Ax, A 2 , Bi and B 2 are invariant when a is replaced 
by 1/a, or when is replaced by 1//3, up to permutations of rows and columns. This means 
that the effective potential (A. 5) is extremal at a — /3 = 1. We only need to check whether 
it is stable there. For 7, there is no such symmetry and the stabilization may occur at a 
nonzero value of 7. A numerical computation has shown that indeed the effective potential 
is stabilized at a — /3 = 1 and at some nonzero 7 as long as m z < m. 



Appendix B Semi-local vortex 

In this appendix, we will review basics aspects of a semi-local vortex (see [51] for a review). 
Usually, a vortex solution is studied by the manifold of minima of a potential V. We call the 
manifold M. mu 



is 



which can be written as 



Mr 



Clonal X G 



global 



local 



x H, 



global 



For example, in the original Abrikosov-Nielsen-Olesen (ANO) vortex, -M m i n = U(l). The 
theory includes only a U(l) gauge group, which parametrizes the minima of the potential. 
When 7r 1 (A / l m i n ) is nonzero there should be either a finite energy vortex or an infinite energy 
one (global vortex) . So we have to check whether or not the kinetic term gives a finite contri- 
bution to the energy. However, nonzero 7Ti (M m m ) does not necessarily imply the presence of 



18 One naturally expects that quantum corrections lift the minimum of the potential unless it is not protected 
by symmetry, so this minimum of the potential generally has the structure G/H, where G and H are the 
symmetry of the theory and that preserved at a vacuum, respectively. 
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a finite energy vortex solution when the theory has some global symmetry and vice versa. A 
vortex generated by symmetry breaking including both global and local symmetry is called a 
semi-local vortex in a broad sense. It includes two types of finite energy vortices: 

• The first type is the same as above: 7Ti(A1 mm ) is nonzero and the symmetry which 
supports the vortex is a gauge group. Therefore the kinetic term does not diverge. Con- 
ventionally, this type of vortex is called the Abrikosov-Nielsen-Olesen vortex although 
this is a semi-local vortex in a broad sense. 

• According to Achucarro, Vachaspati [55], even when 7i"i(.M mm ) = 0, it is possible to 
create a vortex with finite energy. This type can happen only when both local and 
global symmetries exist in a theory. This vortex is called a semi-local vortex. Following 
the convention, we call this type of vortex semi-local. 

In any case, to have a finite energy solution, the kinetic term of a field should approach zero 
faster than 1/r. This condition implies the scalar fields must take values on the same gauge 
orbit at infinity, and the map from infinity to this gauge orbit defines the winding number 
7ri(Ga oca i/ ^locai) • The finite energy solutions including both ANO and semi-local vortices are 
classified by this 7ri(Ga cai/#iocai)- 

Here we will review the explicit construction of the second type vortex by following the 
paper [56] and show a connection to 7r 2 (.M m oduii) in the sigma model limit |57J. A'lmoduii is 
called the vacuum moduli and is defined by 

< / A'lmrn i T Clocal 

M mo duii = j , where L orbit = — . 

-^orbit local 

where L orhit is the gauge orbit in -M mm . The model is a U(l) gauge theory with a Lagrangian 

C = -i(£>„$)t£>"$ - ^($t$ _ rff _ l -F^ , (A.13) 

where $ is an 577(2) doublet and the Lagrangian is invariant under the SU(2) global symme- 
try. The topology of the vacuum manifold is S 3 . Since ^(S 3 ) = 0, one naively expects that 
there is no stable vortex solution. However, one can construct a stable vortex. Since we are 
interested in a static solution, the energy density is given by 

£= h Di ^ + ^^-r) 2 ) 2 + Bf . (A.14) 
2 o 

To keep the energy finite, the energy density at infinity must be equal to zero. This implies 

• The scalar takes values in the minimum of the potential, S 3 . 
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• The first term in the energy density also vanishes. This confines the vev of the scalar 
fields in the gauge orbit S 1 of the local U(l) symmetry in S 3 . 

• The gauge field must be a pure gauge and should be taken so that Dj$ vanishes. 

Therefore, if one identify all the points in the same gauge orbit, the scalar fields take the 
same value in any direction at spatial infinity. Namely we regard as a point. Roughly, 
this means the two-dimensional space transverse to the vortex is identified as S 2 . Now we can 
consider a map from this spatial S 2 to space S 3 /U (1) = S 2 , or more generally, Ggiobai/i^giobai = 
S*?7(2) global /C/(l)diag = S 2 . This is classified by vr 2 (S'f/(2) g i obal /L r (l)diag)- If a configuration is 
trivial under this ir 2 , then we can take a scalar value in spatial CP 1 . Thus, without causing a 



divergence of the kinetic term one can change the configuration to a trivial one 1 ^ It is worth 
noting the correspondence between 7ri(Gi OC ai/-Hiocai) an d ^(-Mmodiuii)- The semi-local vortex 
considered in Gibbons et al. is a one-parameter generalization of the "usual" vortex solution. 
By "usual", we mean that the scalar fields take values in the same gauge orbit everywhere 
on R 2 transverse to the vortex. In this sense, the usual vortex is some limit of a semi-local 
vortex if there exists a one-parameter generalization. On the other hand, semilocal vortices 
can be understood as lump solutions in non-linear sigma models obtained in the low-energy 
limit and whose target space is the vacuum moduli. These lump solutions can be classified 
by ^(Cr global /-^global), so it should agree with 7ri(G ? i OC ai/-ffiocai) — they both classify semi-local 
vortices in different limits. Abrikosov-Nielsen-Olesen vortex is an exception, because it does 
not have any one-parameter generalization. 



Appendix C D-brane and cosmic strings 

In [15], D-brane configurations are presented for all supersymmetry breaking vacua in the 
model studied in section 2 in the main text. In this appendix we would like to identify brane 
configurations for ANO and semi-local vortices. We start with a review of a magnetic brane 
configuration shown in [15] and its slight modification for our purpose. Consider Type IIA 
superstring theory in the flat 10- dimensional Minkowski spacetime with coordinates x ''"' 9 . 
Introduce three NS5 branes. The first one, which we call NSx, is located at x 7,8,9 = and 
extended in the x ''"' 3 and x 4,5 directions (z = x 4 + ix 5 in complex coordinates.). The other 
two are going in directions w = x 7 + ix 8 and x 4,5,9 = but placed at different points in x 6 . We 
call them NS2, NS3. Also we introduce a D6 brane extending in x 4 ' 5 ' 9 directions and placed at 
x 7 ' 8 = 0. For lower energy configurations, we then suspend (Nf — N c — n) D4 branes between 

19 Note that if we weakly gauged S , J7(2) g i b a i, then the configuration is not stable anymore. 
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D6 and NS52, Nq — n D4 branes between NS2 and NS3, and n D4 branes between D6 and 
NS3. The total intersecting brane configuration is shown in Figure [3j 

Before identifying a vortex configuration in this D-brane system, let us comment on the 
gaugino mass and tachyonic direction from the point of view of D-branes. Suppose we embed 
the SM group into SU(Nf — N c — n). This open string mode connecting between Nf — N c — n 
and N c — n corresponds to a messenger. When the distance between these D4 brane becomes 
very close, a tachyonic mode is developed and it causes a reconnection of D4 branes. As was 
shown in the section 2, in this case, the leading order of gaugino masses are non-zero. 

On the other hand, if the SM group is embedded into SU(n), the messenger, which is 
open string connecting two n D4 and N c — n D4 branes, can never be tachyonic. In this case, 
the leading order gaugino masses are vanishing. This fact is consistent with argument shown 
in [23]. 

In a supersymmetric configuration, ANO and semi-local vortices were already studied by 
Hanany and Tong [5B] (see [H] for review). They claimed that a soliton with co-dimension 
two is k D2 branes which are suspended between a NS five brane and D4 branes. The world- 
volume of the D2 branes is extending in x ' 1 ' 2 . We can identify two kinds of such D2-branes. 
One is suspended vertically connected to NSi brane. Since extra D4 branes are not connected 
to the NSi brane, there is no extra degree of freedom. On the other hand, D2 branes going in 
the horizontal direction can have an extra-degree of freedom since Nf — N c — n D4 branes are 
connected to the NS2 brane. We identify this vertex as a semi-local vortex and the other D2 
branes are ANO vortex which we found in the main text. The tension which is proportional 
to the length of the D2 branes is consistent with the symmetry breaking scale of the group 
under which the vortex is generated. 

Finally, we discuss the decay of semi-local vortices. As reviewed in Appendix B, a semi- 
local vortex becomes unstable when one gauges global symmetries. In the brane configuration, 
gauging of non-abelian groups in metastable vacua corresponds to the replacement of the D6 
brane with a NS brane. Once we replace it, a gauge field on D4 brane is introduced, so one 
may say that the end point of D2 can be moved freely and can shrink to zero size. On the 
other hand, if there is no gauge field on D4, moving the end of D2 branes costs infinite energy. 
So it cannot move and the vortex cannot shrink. It would be interesting to explore further 
on this decay process and study the divergence of the kinetic term from the D-brane point of 
view. 
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NS2 

Figure 3: Two kinds of vortices. One is Abrikosov-Nielsen-Olesen vortex which has tension 
fix ~ m 2 - It is D2 brane extending vertical direction. On the other hand, horizontal D2 brane 
is semi-local vortex. The tension is fi? ~ yU 2 . For precise argument one of the NS branes in 
the vertical direction should be replaced D6 for semi-local vortex. 

References 

[1] M. Dine, W. Fischler and M. Srednicki, Nucl. Phys. B 189, 575 (1981); S. Dimopoulos 
and S. Raby, Nucl. Phys. B 192, 353 (1981); 

[2] M. Dine and W. Fischler, Phys. Lett. B 110, 227 (1982); Nucl. Phys. B 204, 346 (1982); 
C. R. Nappi and B. A. Ovrut, Phys. Lett. B 113, 175 (1982); L. Alvarez-Gaume, M. Claud- 
son and M. B. Wise, Nucl. Phys. B 207, 96 (1982); 

[3] S. Dimopoulos and S. Raby, Nucl. Phys. B 219, 479 (1983). 

[4] I. Affleck, M. Dine and N. Seiberg, Nucl. Phys. B 256, 557 (1985). 

[5] M. Dine and A. E. Nelson, Phys. Rev. D 48, 1277 (1993) | arXiv:hep-ph/9303230] . 

[6] M. Dine, A. E. Nelson and Y. Shirman, Phys. Rev. D 51, 1362 (1995) |arXiv:hep- 
|ph/9408384] . 

[7] M. Dine, A. E. Nelson, Y. Nir and Y. Shirman, Phys. Rev. D 53, 2658 (1996) |arXiv:hep- 
|ph/95073"78] . 

[8] M. Dine and J. D. Mason, Phys. Rev. D 78, 055013 (2008) [arXiv:0712.i355] [hep-ph]]. 
[9] K. A. Intriligator, N. Seiberg and D. Shih, JHEP 0604, 021 (2006) [arXiv:hep-th/0602239] . 
[10] K. I. Izawa and T. Yanagida, Prog. Theor. Phys. 95, 829 (1996) |arXiv:hep-th/9602180] . 



26 



[11 

[12 
[13 

[14 
[15 
[16 

[it; 

[18 

[19 
[20 
[21 
[22 

[23 

[24; 

[25 
[26 

[27 



J. de Boer, K. Hori, H. Ooguri and Y. Oz, Nucl. Phys. B 522, 20 (1998) [arXiv:hep- 
|th/9801060] . 



K. A. Intriligator and S. D. Thomas, Nucl. Phys. B 473, 121 (1996) |arXiv:hep- 
th/9603158] . 



H. Murayama, Phys. Rev. Lett. 79, 18 (1997) |arXiv:hep-ph/9705271] . 

M. Dine and J. Mason, Phys. Rev. D 77, 016005 (2008) |arXiv:hep-ph/0611312| ; M. Dine 
and J. D. Mason, |arXiv:0712.1355| [hep-ph], 

M. Ibe and R. Kitano, Phys. Rev. D 75, 055003 (2007) | arXiv:hep-ph/0611111 |; M. Ibe 



and R. Kitano, JHEP 0708, 016 (2007) |arXiv:070 5.3686 [hep-ph]]. 
R. Kitano, H. Ooguri and Y. Ookouchi, Phys. Rev. D 75, 045022 (2007) larXiv:hep- 



ph/0612139] . 



C. Csaki, Y. Shirman and J. Terning, JHEP 0705, 099 (2007) |arXiv:hep-ph/0612241] . 
R. Essig, J. F. Fortin, K. Sinha, G. Torroba and M. J. Strassler, JHEP 0903, 043 (2009) 



|arXiv:0812.3 213 [hep-th]] 



H. Ooguri and Y. Ookouchi, Phys. Lett. B 641, 323 (2006) |arXiv:hep-th/0607183|. 



S. Franco, I. Garcia-Etxebarria and A. M. Uranga, arXiv:hep-th/0607218 



I. Bena, E. Gorbatov, S. Hellerman, N. Seiberg and D. Shih, arXiv : hep-th /0608157[ 



K. A. Intriligator and N. Seiberg, Class. Quant. Grav. 24, S741 (2007) [arXiv:hep- 



ph/0702069] 



R. Kitano, H. Ooguri and Y. Ookouchi, [arXiv: 1001 .4535 [hep-th]. 



K. I. Izawa, Y. Nomura, K. Tobe and T. Yanagida, Phys. Rev. D 56, 2886 (1997) 
|arXiv:hep-ph/9705228| . 



Z. Komargodski and D. Shih, JHEP 0904, 093 (2009) (arXiv:0902.0030 [hep-th]]. 



A. Giveon, A. Katz and Z. Komargodski, JHEP 0907, 099 (2009) |arXiv:0905.3~3~8~7 [hep- 
th]]. 



M. Ibe, Y. Nakayama and T. T. Yanagida, Phys. Lett. B 649, 292 (2007) larXiv:hep- 
ph/0703110] . 



27 



[28] M. Ibe, K. I. Izawa and Y. Nakai, Phys. Rev. D 80, 035002 (2009) |arXiv:0812.4089 



[hep-ph]]; M. Ibe, K. I. Izawa and Y. Nakai, |arXiv: 0907.2970 ) [hep-ph]. 

[29] 



[30 
[31 
[32 

[33 
[34; 
[35 
[36 
[37 

[38 

[39 

[40 

[41 
[42 

[43 
[44; 
[45 



S. A. Abel, J. Jaeckel and V. V. Khoze, arXiv:0907.0658 [hep-ph]. 



K. Hanaki and Y. Ookouchi, arXiv: 1003.5663 [hep-ph]. 



Y. Nakai and Y. Ookouchi, arXiv: 1010.5540 [hep-th]. 



D. Green, A. Katz and Z. Komargodski, Phys. Rev. Lett. 106, 061801 (2011) 
|arXiv:1008.22l5l [hep-th]]. 



R. Auzzi and A. Giveon, JHEP 1010, 088 (2010) |arXiv: 1009.1714] [hep-ph]]. 
M. Sudano, |arXiv:1009.2"086| [hep-ph]. 



T. Banks, N. Seiberg, Phys. Rev. D83, 084019 (2011). |arXiv:1011.5120 [hep-th]]. 



T. Banks and H. E. Haber, JHEP 0911, 097 (2009) |arXiv: 0908.2004 [hep-ph]]. 

A. D. Dolgov and L. A. Ponomarev, Yad. Fiz. 5, 164 (1967) [Sov. J. Nucl. Phys. 5, 114 
(1967)]. 

A. D. Dolgov, Yad. Fiz. 7, 394 (1968) [Sov. J. Nucl. Phys. 7, 255 (1968)]; C. L. Basham 
and P. K. Kabir, Phys. Rev. D 15, 3388 (1977). 

M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 (1986); For a recent review, 
W. Buchmuller, R. D. Peccei and T. Yanagida, Ann. Rev. Nucl. Part. Sci. 55, 311 (2005). 



M. Eto, K. Hashimoto and S. Terashima, JHEP 0703, 061 (2007) |arXiv:hep- 
|th/0610042] . 



D. Tong, arXiv:hep-th/0509216 



E. J. Copeland and T. W. B. Kibble, Proc. Roy. Soc. Lond. A 466, 623 (2010) 
[arXiv:0911.1345l [hep-th]]. ; T. W. B. Kibble, [arXTv:astro-ph/0410073[ 



W. H. Zurek, Nature 317, 505 (1985). 

H. Murayama and J. Shu, Phys. Lett. B 686, 162 (2010) [arXiv:0905.1720| [hep-ph]]. 



A. Vilenkin and E. P. S. Shellard, Cambridge University Press (1994). 



28 



[46] M. Wyman, L. Pogosian and I. Wasserman, Phys. Rev. D 72, 023513 (2005) [Erratum- 
ibid. D 73, 089905 (2006)] |arXiv:astro-ph/0503364] . 

[47] A. Vilenkin, Phys. Rev. D 23, 852 (1981). 



[48; 

[49 

[50 

[51 

[52 
[53 
[54 

[55 
[56 

[57 
[58 



V. Vanchurin, K. D. Olum and A. Vilenkin, Phys. Rev. D 74, 063527 (2006) |arXiv:gr- 



qc/0511159]. 



C. Ringeval, M. Sakellariadou and F. Bouchet, JCAP 0702, 023 (2007) [ arXiv:astro- 
ph/0511646] . 



C. J. A. Martins and E. P. S. Shellard, Phys. Rev. D 73, 043515 (2006) | arXiv:astro- 
ph/0511792] . 



M. Kawasaki, K. Miyamoto and K. Nakayama, Phys. Rev. D 81, 103523 (2010) 
[arXiv:1002.0652l [astro-ph.CO]]. 



M. R. DePies and C. J. Hogan, Phys. Rev. D 75, 125006 (2007) |arXiv:astro-ph/0702335] . 



T. C. Collaboration, jarXiv: 11 02.2181 [astro-ph.CO]]. 

A. Achucarro and T. Vachaspati, Phys. Rept. 327, 347 (2000) [Phys. Rept. 327, 427 
(2000)] |arXiv:hep-ph/9904229| . 

T. Vachaspati and A. Achucarro, Phys. Rev. D 44, 3067 (1991). 

G. W. Gibbons, M. E. Ortiz, F. Ruiz Ruiz and T. M. Samols, Nucl. Phys. B 385, 127 
(1992) |arXiv:hep-th/9203023] . 

M. Hindmarsh, Nucl. Phys. B 392, 461 (1993) |arXiv:hep-ph/9206229] . 
A. Hanany and D. Tong, JHEP 0307, 037 (2003) |arXiv:hep-th/0306150| . 



29 



